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SUMMARY 

A program for the computation of neutron ages from experimental activation data is 
described. The quadrature schemes and error analyses are detailed along with some 
curve-fitting techniques. A ’’best-fit” derived from statistical decision theory is de- 
scribed. 



INTRODUCTION 

One important parameter found in the critical equation (eq. (1)) of age -diffusion 
theory is the neutron age r (ref. 1): 
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where K is the multiplication factor, B the buckling, and if the diffusion length. The 
neutron age t is related to the mean-square distance traveled by a neutron from its 
point of generation to its point of thermalization (ref. 1) 
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where is the initial energy of the neutron and E is some lower energy generally in 

the thermal region. 

2 

Inasmuch as the value of r is impossible to measure directly, the saturated activ- 
ity A 0 (R), which is directly proportional to the neutron flux distribution (p( E, R) is mea- 
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sured by foil activation methods. The value of r may then be obtained by familiar 
averaging techniques; that is, 

R 4 A g (R)dR 

(3) 

R 2 A g (R)dR 

and the neutron age may be obtained from equation (2). 

The aforementioned set of integrations or averaging process, when done by hand, re- 
sults in a neutron age subject to errors introduced by the individual calculating this age. 

In fact, the error introduced by hand quadrature might be as large as several square 
centimeters. Another problem associated with hand calculation is the assessment of the 
error in the age, that is, the sum of the experimental and calculational errors. 

The ELLI code has been written to reduce the time necessary to calculate a neutron 
age and to provide a definite schema of procedural choices. It thus allows convenient 
comparison of several fitting and error analysis procedures. 

In the analysis, the problem is defined, the numerical methods used are explained, 
and various approaches to the error assessment are presented. 

SYMBOLS 

saturated activity, counts/min 
-2 

buckling, cm 
cofactor array of J | c rg | j 

coefficient array of normal equations of weighted least-squares curve fit 
determinant of array ||c rg || (see eq. (27)) 
breakpoint between region I and exponential region (fig. 1), cm 
lower energy of neutron, eV 
initial energy of neutron, eV 
F test ratio 

polynomial order 

panel width of trapezoidal integration, cm 
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average panel width, cm 
multiplication factor 
order of approximating polynomial 
sum of ARj (see eq. (20)), cm 
evaluated by equation (21) 
number of tabular values 
degree of polynomial 
proportionality constant 

slowing down density, neutrons/cm • sec 
distance from source, cm 
R. corresponding to minimum value of 

o 

mean-square distance traveled by a neutron, cm 
estimate of a 2 
weighting function of y(R) 
th. 

i value of approximating polynomial 
index number of breakpoint 
residual of R. 
error estimate 

macroscopic absorption cross section of indium, cm 

variance of curve fit 
neutron age, cm 2 

o 

neutron flux distribution, neutrons/cm 


PROBLEM ANALYSIS 

The Fermi age is related to the second moment of the slowing down density q(R, E) 
by equation (2). Experimentally, however, the flux age is measured by finding the spa- 
cial distribution of neutrons at a specific energy. For example, the activity, if satura- 
tion is assumed, of a cadmium-covered indium foil at a distance R from the source is 
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Cd cutoff 


Since the value of the activation of the indium foil is nearly zero at all energies ex- 
cept that of the indium resonance, 1. 44 electron volts, the flux age for the indium reso- 
nance neutrons, can be expressed as 


T ( E i> E i 44 e y) 


r 2 < E i' E 1.44eV> 
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R 4 A Q (R)dR 
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R 2 A g (R)dR 
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in spherical coordinates. This equation is then evaluated by the ELLI code. The methods 
by which these integrals are evaluated are outlined in the Numerical Analysis section. 

The calculation of the errors associated with these methods is shown in the Error Analy- 
sis section. 


Numerical Analysis 


For the purpose of integration, the experimentally determined activation times dis- 
tance squared curve shown in figure 1 is separated into two portions. The shape of 
curve I has not been analytically determined, whereas curve n is assumed to be exponen- 

ential (ref. 2). The breakpoint or division 
between curves I and n may be determined 
by the program or the experimenter. E the 
program option is selected, the derivatives 
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of A„(R)R with respect to R 
s 

lated by a difference method: 




A s< R l + l> R ll - A s < R i> R i 
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l i 


Figure 1. - Representative curve of distance squared times -saturated 
activity against distance, with division point shown. 
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'The value of R. corresponding to the minimum value of the differences, that is, the 
largest negative value of D^, is chosen as the breakpoint between curve I and the expo- 
nential region. H this specific value R g is not suitable, the breakpoint value can be 
chosen by the experimenter as in opticm 4 (see appendix). 

Curve I is then integrated in any of several fashions that are left as options in the 
code. The first option available for the integration of curve I is the trapezoidal integra- 
tion. This option should be used if many finely spaced points are to the left of the division 
point Rg. The symbolic rule used for integration is 
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The contributions to the integrals are printed out, point by point, for this option. Sample 
computer output and a listing of the program are shown in the appendix. 

The second option available is a least-squares polynomial curve fitted to the data. 
This option also allows the specification of the significance of the fit; that is, a certain 
confidence level for the polynomial fit may be inserted as pregram input. The program 
first attempts a quadratic fit to the curve I data, and then, if the confidence level crite- 
rion is not fulfilled, the program uses a cubic fit. This process is repeated, with the de- 
gree of the polynomial raised after each step, until the significance criterion is satisfied, 
or until one of two other terminating conditions is found to exist. The two other termi- 
nating conditions are when (1) a tenth-degree polynomial is reached or when (2) the de- 
gree of the polynomial reaches one less than the number of points to be fitted. When 
either terminating condition is obtained, the appropriate action is taken. S condition (1) 
exists, the tenth-degree polynomial is integrated analytically; if condition (2) exists, the 
(n-l) th polynomial is likewise integrated. The test of significance is computed by a 

method following that of reference 3. First s 2 , the estimate of the variance for the 
2 

curve fit a , is computed: 
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If the estimate of the variance computed for the g order polynomial by c r is denoted, 

th © 

then the test for significance at the p confidence level, or probability of rejection, is 
found to be 


( n - g) 6 e_i - (n - g - 1 )o 2 
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where F/„ \ is the "well-known” F ratio of statistical analysis. If the change in 

v 2' , 

variance is not significant at the specified level, the g tft polynomial is accepted and inte- 
grated by analytic means over the range from zero to Rg. That is, if the probability 
level chosen is 0. 500 and the probability that a cubic fit is better than a parabolic is only 
0.400, the parabolic is accepted and integrated analytically. 

The third option available is an averaged parabolic fit similar to that of reference 4. 
This option fits one parabola to the first three and another to the last three of four suc- 
cessive points. The coefficients thus generated are averaged and the curve is integrated 
in the overlapping region as shown in figure 2. Li this option, the integration is done in 
small portions; for example, in figure 2 the integration would be from Rg to R^: 



This option produces a considerably better approximation to the integral than with the 
trapezoidal option inasmuch as curvature is allowed between points. 

The exponential region of the curve, curve n, is fitted by a weighted least-squares 

O 

exponential approximation. The A (R.)Rr values are transformed logarithmically by 

Sll 
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A S (RIR 2 , (counts)! sq cm)/min 


Points 

used 




Normalized weights are assigned proportional 
to the inverse of the square of the standard devia- 
tion or variance as determined by the counting sta- 
tistics of a data point. Then a weighted least- 
squares linear fit is made to obtain the coefficients 
Sq and a^, such that 

y i ^Y i = a x x. + a Q (13) 

and subject to the condition that the weighted sum 
of the residuals squared S is minimized: 

N N 

s = E W A> 2 s £ w i (y i - y/ (14) 
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Finally, the variables Sq and a^ are transformed back from the semilogarithmic coor- 
dinates by the transformation 



(15) 


The approximation 
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of curve n is integrated analytically from R g to 00 . The contribution to the integrals 
a and 0 appear in the output from the program. Option 4 allows the value of 6 to be 
specified as part of the program input. 
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Error Analysis 


The remaining options of the ELLI code provide for error estimation of the various 
integration options. 

The error estimate of a trapezoidal integration is given by (ref. 5) 

= (17) 

where h is the panel width of the trapezoidal integration. Option 5 of ELLI uses h, the 
average panel width, instead of h and the derivatives of a parabolic fit to the data in the 
range (0, R g ) for the value of d/dR|A s (Rg)RgJ. The fractional errors 6^ and of the 
two integrals (eqs. (7) and (8)) are then combined over the range (0,R fi ) while indepen- 

V ~2 2 

9 1 + #2 is assumed. This allows a fractional error to be as- 
signed to the integration of curve I by the trapezoidal rule. 

Option 6 assigns an error to the integration of option 2, which is an integration of a 
g** 1 degree polynomial fit. The error bound may be expressed as (ref. 6) 


where 
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in (0, R 6 ) 


l = r 6 -° 


(18) 


(19) 


and W(R) is the weighting function of y (R). As the weighting function is relatively con- 
stant over curve I, a value of 1 shall be assigned to W(R.), and / 6 W(R)dR shall be 
i=6 1 Jo 

replaced by ^ AR^ = L. Upon substitution, the following is obtained: 
i=l 


2 g (g + 1)! 


( 20 ) 
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" In order to obtain an approximation of M, a (g + 2)th degree polynomial is fitted to the 
data, and the (g + l)th derivative is taken and evaluated at R g : 


M « 


d ( g +1 ) 

djjfe +1 ) 


y(R fi ) 


( 21 ) 


For both integrals is then computed, and a fractional error is assigned to the integra 
tion of curve I in the same m ann er as in option 5. 

Option 7 is again an error ass ignm ent - in this case, to the option 3 integration of 
curve I. The error is given (ref. 5, p. 385) by 


*3 




( 22 ) 


because option 7 is similar to an unequally spaced Simpson’s integration. Again, h is 
approximated by h, A g (R)R is approximated over the interval (0, Rg) by a fourth degree 
polynomial y(R), and the third derivative of this approximating polynomial is evaluated 
at Rg. The fractional error associated with the two integrations over curve I is printed 
as in the previous options. 

Finally, the error associated with the integration of the exponential portion, curve H, 
is calculated. This error calculation is described separately because it is computed in 
the same manner in options 5 to 7. 

Reference 6 (pp. 261 to 269) shows that the error associated with the coefficients of 
a weighted polynomial approximation to tabular data may be represented by 



For the case in question (see eq. (14)), the error formulae reduce to 
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where 


Q = 


:[w(R.)Tj(R i ) 2 ] 
N - 1 


D = [sW(R i )] [SW(R.)R 2 ] - [SWCR^R.] 2 (26) 

and the summation runs from 6 to the number of data points. Hence, if the transforma- 

a 0 

tion /3 = a^ and a = e is used, to gain the exponential approximation, the errors 
transformas 6/3 = 6a^ and 5a» aSa^. These errors are then propagated through the 
analytic integrations as standard errors. 


CONCLUDING REMARKS 

The ELLI code has been operative at Lewis for approximately 1 year. Experience 
gained by the operation of the program throughout this period has shown that all options 
are used about equally. The particular option chosen is to some degree dependent on the 
data and its quality. In general, the ages obtained through the use of the first three op- 
tions vary by only a few tenths of a square centimeter. If the variance of these ages is 
much larger, more data are probably necessary. 

Operation times for this code executed on the Lewis 7094-H computer are character- 
istically of the order of hundredths of a minute. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 10, 1966. 
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APPENDIX - INPUT TO ELLI 


CARD 1 FORMAT (A6) 

bbbbbb in first six columns (if not blank a program description is output) 

CARD 2 FORMAT (711) 

Card columns 1-7 

CARD 3 FORMAT (A6, 4X, A6, 4X, 8A6) 

Card columns 1-6 ASbbbb H activations are to be read. 

o 

ASR** 2 If activations times R are to be read. 

Card columns 11-16 bbbbbb H data are in centimeters. 

INCHbb E data are in inches. 

Card columns 21-68 Label for experiment. 

CARD 4 FORMAT (18, 12X, F10. 3) 

Card columns 1-8 Number of data points <100. 

Card columns 21-30 Probability level of significance, 0. 500 s p 

<0. 995; needed only if option 2 is specified. 

CARD 5 FORMAT (3F 15. 4) 

Card columns 1-15 Rj in inches or centimeters according to CARD 3. 
Card columns 16-30 A g (Rj) or A g (R.)R? according to CARD 3. 

Card columns 31-45 ajA. g (R.)] or o|A g (Rj)R?J according to CARD 3. 

LAST CARD FORMAT (18) 

B option 4 is chosen. 

Card columns 1-8 6 is index of division point Rg. 

For more than one set of data, repeat starting with CARD 2 of INPUT TO ELLI. 


0 in column n if n th option to be skipped; 1 if n th option to 
be used. Options 5, 6, and 7 may be used only if the corre- 
sponding option 1, 2, or 3 has been used. 
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ELLI LISTING 


EOUIVALENCEIAZERO, ALPHA) , ( T( l ) «X( 1) ) ,IT!101),Y!l)) ,!T(201),STD(1)) 
l,!T1401),DEV!l)) 

DI MENS IONOE V< 100) , STD! 100) ,T( 100,5) ,X( 100) ,XE( 100) , Y ( 100 ) , LOGTR! 10 
10 ) «R21 100) • DERL ( 100), C 11 100),C2!100) ,C31100) ,BLEKi ( 8) , LABE ( 8 ) , A1 ( 2 
20) ,A2( 20),A(20),ARB1(20) ,ARB2( 20) ,Y2! 100) ,XTEMP! 100) , TEMPI 100), WE I 
3GH! 100) 

REAL LOGTR 

REAL INTR2, INTR4, INTRP2, INTRP4 

INTEGER OPT1,OPT2,OPT3,OPT4,OPT5,OPT6,OPT7 

INTEGER EXPON, ZERO,OATAK , ASRSQ,AS ,ORO,ORON,OROP1 ,0R0P2 

COMMON IVOZEO 

1000 REAOt 5 , 10 10 ) HOW,BLEK,BLEKi 

1010 FORMAT ( A6, 4X ,A6, 4X , 8A6 ) 

1011 FORMAT ( 71 1 ) 

10 WRITE (6,20) 

20 FORMAT! 1H1,45X,36H AGE MEASUREMENT BY FOIL ACTIVATION. //) 

1012 READ! 5, 10 l 1 )OPT l, 0PT2,0PT3,QP T4,0PT5 ,0PT6 ,0PT7 

22 00 25 JB = 1,100 

23 00 24 KB = 1,5 

24 TUB, KB) = 0.0 
LOGTR(JB) = 0.0 
R2( J B) = 0.0 
OERL(JB) = 0.0 
Ci(JB) = 0.0 
C21JB) - 0.0 
C3( JB) = 0.0 
Y2( JB) = 0.0 
XTEMP(JB) » 0.0 
TEMP(JB) = 0.0 

25 WEIGHIJB) = 0.0 

26 DO 27 JK = 1,20 
A1UK) = 0.0 
A2(JK> = 0.0 

A( JK ) = 0.0 

AR8KJK) = 0.0 

27 ARB2( JK ) = 0.0 
35 KLUG = 0 

1015 DATA BLAK.ASRSQ, AS/0606060606060, 0216251545402, 0216260606060/ 

1016 OATA INCH/0314523306060/, ZERO/O/ 

1020 IF1H0W-BLAK)1030, 1050,1030 

1030 CALL HOWAGE 
1040 STOP 

1050 READ( 5,1010) DATAK,LUM,LABE 
1060 IF ( DAT AK -A SR SO) 1070, 50, 1070 
1070 IF(0ATAK-AS)1080, 1110, 1080 
1080 WRITE (6,1090) 

1090 FORMAT! 20X.67H DATA FORMAT CARD INCORRECT OR MISPLACED. SHOULD BE 
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i ASR**2 OR AS. ) 

1100 STOP 
1110 KLUG a 1 
1120 GO TO 50 

1130 DO 1140 IPER = 1* IDATA 

1135 STDt IPER) = STD I IPER)*XIIPER)*X(I PERJ 

1140 Y( iPER ) = XI IPER)*Xt IPER)*Yf IPER) 

1150 GO TO 90 

50 READt 5*60 ) IDATA, EPSILO.PE 

60 FORMAT ( I8»2X»F10.3,F10.3) 

70 READt 5,80 ) 1 1 T€ IK, JUP) » JUP*1 *3),IK=1,IDATA) 

00 FORMAT C3F15.4I 
9000 IF10PT4T81, 81,9010 
9010 READt 5, 90201 IPOL 
9020 FORMAT (181 

81 I F t LUM— INCH 1 87, 82, 87 

82 DO 84 1 - 1,100 

84 XU J = XI I)*2. 54 

C IF INCHES USED NO ASR**2 CARD 

85 IFl LUM- INC H.EQ .ZERO . AND. DATA K-ASRSQ.EQ. ZERO) GO TO 86 
GO TO 87 

86 NRITEI6, 11601 

1160 FORMAT t 18H ERROR CONDITION 4) 

STOP 

87 I Ft KLUG I 1130,90, 1130 

90 IMAX«1 

9400 DO 9410 I = 1. IDATA 

9410 Y 2 tn * y<n*xtn*xm 

91 IFtXt 11)92, 100,94 

92 WRITEt 6, 93) 

93 FORMAT t 18H ERROR CONDITION 1) 

STOP 

94 IOATA * IDATA ♦ 1 
9030 IFt0PT4)9050, 9050, 9040 
9040 I POL * l POL ♦ l 

9050 I DM2 a IDATA - 2 

DO 95 I = ZERO, 10M2 
IPl *1+1 
IADAP * I DATA- IPl 
I ADAS * IDATA -I 
X( I ADAS ) = Xt IADAP ) 

Yt I ADAS) = Yt IADAP) 

Y2< IADAS) * Y2t IADAP) 

95 STDt IADAS) = STOt IADAP) 

Y2t l) = 0.0 

Xtl) * 0.0 
Yt 1) » 0.0 
STDt l) * STDt 2) 

100 IF(0PT4)9O60, 9060,9070 
9070 IZORK = I POL 
9080 GO TO 2000 
9060 DO 110 IIP = 2, IDATA 

110 DERL (IIP) = (Yt IIP)— YtIIP-l))/tXtII P)— X t I IP— 1 1 ) 

115 DERL ( 1 ) = 0.0 
120 MIN = 1 



130 DO 150 IZORK = 1,IDATA 

140 I F ( DERL ( MI N I. GT .DERL 1 IZORK ) ) MI N = IZORK 
150 CONTINUE 
160 I POL = MIN 
2000 IFI0PT1 12020,2010,2050 
2010 GO TO 3000 
2020 WR I TE I 6, 2030) 

2030 FORMAT! 18H ERROR CONDITION 2) 

2040 STOP 

2050 OPT 1 = OPT l- 1 

2051 WRITE (6,2060) LABE 

2060 FORMAT! 54X,9H OPTION l,iOX,8H LABEL =,8A6//) 

2070 WRITE! 6,2080 1 

2080 FORMAT ( 35X, 2H 2/13X.2H R,17X,4H A R,llX,14H ST. OEVIAT I0N.9X, 12H P 

I ART I AL SUM , 10X, 12H PARTIAL SUM/33X.2H S.34X,14H OF INTEGRAL i,8X, 

2 14H OF INTEGRAL 2) 

2090 SUM 1 = 0.0 
2100 SUM2 = 0.0 

2110 SUM 1 = ( 1.0/3. 0)*X(2)*Y(2) 

2120 SUM2 = ( 1.0/3. 0 )* Y( 2 l*X ( 2 ) **3 
2130 WRITE(6,2140) X! 1 ), Y! 1) ,STD( 1 ) 

2140 FORMAT! 7X.2PE15.7, 4X ,2PE 15. 7, 6X, 2 PE 15. 7) 

2150 WRITE! 6,2210) X! 2 ) , Y I 2 ) , STD ( 2 1 , SUM1 , SUM2 
2160 IFUD ■ I POL - 1 
2170 DO 2200 I * 3, IFUD 

2180 SUM1 = SUM1 ♦ 0.5*(YI I )+Y( I-i) )*(X! I )-X( 1-1) ) 

2190 SUM2 = SUM 2 ♦ 0.5*( ( Y( I 1*X ( I ) **2) ♦ ! Y( 1-1 ) *X< I- l ) **2 )) * ( X ( I )-X ( I- 1 ) 
1 ) 

2200 WRITE(6«2210)X( I ) , Y( I), STD (I ) ,SUMl,SUM2 

2210 F0RMAT!7X,2PE15.7,4X,2PE15. 7, 6X, 2PE15. 7 ,6 X,2PE 15. 7 ,8X, 2PE 15.7 ) 

2220 SUM 1 = SUM 1 +0. 5* ( YUPOL) + Y 1 1 POL-1 1 ) * ! X! I POL 1 -XU POL-l ) ) 

2230 SUM2 = SUM2 ♦ 0.5*11 Y ( IPOL I *X( I POL » **2+Y II POL-1 1 *Xt I POL-l ) **2 ) *( XI 

I I POL I— X ( I POL- 1)1 ) 

2240 WRITE ( 6,2250 ) X ( IPOL ) , Y ( I POL ) , STD ( I POL) , SUMl , SUM2 
2242 TRA1 = SUMl 
2244 TRA2 = SUM2 

2250 FORMAT! 3X , 4H**** , 2PE 15. 7,4X, 2PE15. 7 ,6X ,2PE 15. 7 ,6X, 2PE l 5. 7 . 8X, 2PE15 
1.7, 4H****///) 

2260 ASSIGN 2280 TO ICOM 
2270 GO TO 8000 
2280 SUMl = SUMl + INTRP2 
2290 SUM2 = SUM2 + INTRP4 
2300 AGET = SUM2/I 6. 0*SUM l ) 

2310 WRITE (6,2320) 

2320 FORMAT! 43X.32H EXPONENTIAL REGION CALCULATI 0NS//35X.2H 2/13X, 

12H R,17X,4H A R,11X,14H ST. DE VI ATI0N.8X , 13HAPPR0X IMAT ION, 7X, 

2 16H PERCENTAGE DEV./33X,2H S) 

2330 DO 2360 I LLK * IPOL,IDATA 
2340 APR = ALPHA*EXP(BETA*X( ILLK) ) 

2350 PDEV = ( ( Y ( ILLK )— APR l/STD! ILLK) )*100.0 

2360 WRITE (6,2370) X! ILLK ) ,Y( ILLK) , STD! I LLK) , APR, PDEV 

2370 FORMAT! 7X, 2PE 15. 7,4X , 2PE 1 5. 7 , 6X.2PE 1 5. 7 ,6X, 2PE15. 7 ,7X,2PE15 .7 » 

2380 EXLEN = 1. O/BETA 

2390 WRITE ( 6, 2400 ) AGET , A ZERO,E XLEN, I NTRP2 , I NTRP4 

2400 FORMAT ( 42X , 38H AGE COMPUTED BY TRAPEZOIOAL OPTION IS/45X,2PE15. 7, 
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1 » 1 1HSQUAR E CMS.///20 X, 8H AZERO = , 2PE15. 7 ,20X ,20H RELAXATION LENG 
2TH =,2PE15.7//15X,31H INTEGRAL 1 OVER EXPONENTIAL IS.20X, 

33IH INTEGRAL 2 OVER EXPONENTIAL IS/ 22 X, 2 PE 15. 7 »3 AX ,2 PE 15.7 > 

2410 IFI0PT5) 2420, 2450, 2460 
2420 WRITEI6. 24301 

2430 FORMAT ( 18H ERROR CONDITION 31 

2440 STOP 

2450 GO TO 3000 

2460 0PT5 = 0PT5 - 1 

2470 GO TO 5000 

3000 IF!OPT2)3010, 3040, 3050 

3010 WRITE16, 3020) 

3020 FORMAT! 18H ERROR CONOITION 7) 

3030 STOP 
3040 GO TO 4000 
3050 OPT 2 « 0PT2 - 1 
3060 WRITE(6»3070) LA8E 

3070 FORMAT! 1H1,45X*36H AGE MEASUREMENT BY FOIL ACTIVAT I0N.//34X, 9H OP 
IT ION 2, 10X, 8H LABEL -,8A6//} 

3080 WRITE! 6,3090) 

3090 FORMAT ( 40X, 2H 2/20X.2H R, 15X.4H A R,15X,19H STANDARD DEVIATION, 

1 14X* 14H APPROXIMATION, 7X.IIH PERCENT AGE/38X ,2H S,70X,10H DEVIATION 
2 ) 

3100 ORD = 2 

3110 ORDN = 3 

3120 FORD = FLOAT! ORD) 

3130 FORDN - FLOAT! ORDN) 

3140 NN = I POL 

3150 FNN = FLOAT! IPOL) 

3160 DO 3400 IBIS = 2,9 

3170 CALL POYFZ! X, Y, I PCX. ,ORD,A l ) 

3180 CALL POYFZ!X,Y, IPO. , ORDN, A2) 

3190 SIGMA1 = 0.0 
3200 SIGMA2 * 0.0 
3210 DO 3260 IKIS-l.NN 
3220 YAP =A1! 1) 

3230 YORP = A2!l) 

3240 DO 3260 (CIS - 1,10 

3250 YAP * YAP ♦ All IC IS+1 )*X( IKI S I**ICI S 
3260 YORP * YORP + A2{ ICI S ♦ 1 )*X( I KI S)**ICI S 
3270 SIGMA1 » SIGMA1 ♦ IYIIXIS) -YAP)**2 
3280 SIGMA2 = SIGMA2 ♦ (Y! IKIS)— Y0RP)**2 
3290 SIGMA1 = SIGMAl/iFNN-FORO-l.O) 

3300 SIGMA2 * SIGMA2/! FNN-FORDN-l.O) 

3305 IF! SIGMA1.LE.O.O.OR.SIGMA2.LE.O.OJGO TO 3410 
3310 NEST = NN -ORD 

3320 IF! !U FNN— FORDN ) *SIGMA1- ( FNN— FORDN- 1 . ) +SIGMA2 ) /SIGMA2 ).LT.FTES!1,N 
1EST.PE) ) GO TO 3410 
3330 IF!NN-0RD-3)3340, 3340, 3370 
3340 WRITE! 6,3350) 

3350 FORMAT! 24H TERMINAL COND. 1 EXISTS) 

3360 GO TO 3410 
3370 ORD = ORD +1 
3380 ORDN » ORDN * 1 
3390 FORD * FLOAT I ORD) 

3400 FORDN = FLOAT! ORDN) 
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3410 XI = X! IPOL l 
3420 XTEMP = 1.0 
3430 XTEH2 = XI * XI 
3440 INTR2 = 0.0 
3450 INTR4 = 0.0 
3455 0RDP1 = ORD + 1 
3460 DO 3500 I = 1.0RDP1 
3470 XTEMP = XTEMP*>X1 
3480 XTEM2 = XTEM2 * XI 

3490 INTR2 = INTR2 + XTEMP*A1! I l/FLOAT ! I ) 

3500 INTR4 = INTR4 + XTEM2 * A 1 ( I ) /FLOAT 11+2) 

3510 APPRO = 0.0 
3520 PDEV = 0.0 
3530 00 3580 I = i.IPOL 
3540 00 3560 J = 1, ORDPi 
3545 JMi = J-l 
3548 IF! JM1)3550, 3550,3555 
3550 XJTJM1 » 1.0 
GO TO 3560 

3555 XJTJM1 = X( l)**JMl 

3560 APPRO = APPRO «• A1(J)*XJTJH1 

3570 PDEV = UYII) - APPRO »/STDC I )) *100. 0 

3575 WRITE (6, 3590 )X( I), Y{ I l,STD( I ), APPRO, PDEV 

3580 APPRO = 0.0 

3590 FORMAT ! 13X, El 5.7, 5X, E 15.7 , 10X, E 15. 7 ,15X, E15. 7 , 5X, 2PE 15.7 ) 

3600 WRITE16, 3610) 

3610 FORMAT 4 1H , 1H*,4X, 1H* , 4X , 1H*, 4X ,1H* ,4X , 1H* ,4X , 1H* , 4X , 1H* ,4X , 1H*, 
14X, 1H*» 4X , 1H*,4X.1H*,4X,1H*,4X,1H*,4X,1H*,4X, 1H*,4X,1H*,4X,1H*, 
24X, 1H*,4X , 1H*,4X, 1H*,4X,1H*,4X,1H*,4X,1H*,4X,1H*,4X ,IH*,4X,1H*) 

3614 WRITE! 6,3615) SIGMA1,PE 

3615 FORMAT! 10X, 17H SIGMA1 SQUARED =,E15. 7 .20X.24H LEVEL OF SIGNIFICANC 
IE =* E 10. 4 ) 

3620 ASSIGN 3640 TO ICOM 
3630 GO TO 8000 
3640 WRITE! 6,2320) 

3650 HSU1 = INTR2 ♦ INTRP2 
3660 HSU2 = INTR4 + INTRP4 
3670 HAGE * HSU2/1 6.0*HSU 1 ) 

3680 DO 3710 ILLK* IPOL , I DATA 
3690 APR = ALPHA*EXP I BETA*X! ILLK) ) 

3700 PDEV =! !Y! ILLK)— APR) /STD! ILLK) 1*100.0 

3710 WR I TE! 6,2370 ) XI ILL K ) , Y ! ILLK ) , STD! I LLK) .APR, PDEV 

3720 EXLEN = 1.0/8ETA 

3730 WRITE! 6,3740) HAGE , A ZERO.E XLEN, I NTRP2 , I NTRP4 , ORO 

3740 FORMAT! 42X.38H AGE COMPUTED BY POLYNOMIAL OPTION IS //45X.2PE 15.7, 
1,1 1HSQUARE CMS.///20X, 8H AZERO =, 2PE15.7 ♦ 20X ,20H RELAXATION LENG 
2TH =, 2PE15.7//15X.31H INTEGRAL 1 OVER EXPONENTIAL IS.20X, 

331H INTEGRAL 2 OVER EXPONENTIAL I S/22X.2PE1 5. 7 ,34X ,2PE 15 .7/ 
450X.21H DEGREE OF EQUATION =,I3) 

3750 IF(QPT6)37 60, 3790, 3800 
3760 WR I TE(6,3770) 

3770 FORMAT! 18H ERROR CONDITION 9) 

3780 STOP 
3790 GO TO 4000 
3800 0PT6 * OPT 6 - 1 
3810 GO TO 6000 
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' . 4000 I F ( OPT 3 ) 40 i 0* 4040* 4050 
4010 WRITE(6. 40201 

4020 F0RHATH8H ERROR CONDITION 51 

4030 STOP 

4040 GO TO 1012 

4050 0PT3 = 0PT3 - 1 

4060 WRITEI6, 4070) LABE 

4070 FORMAT ( 1H1.45X.36H AGE MEASUREMENT BY FOIL ACTI VATI0N.//34X, 9H OP 
IT ION 3* 10X» 8H LABEL *.8A6//) 

4080 WRITE! 6, 2080) 

4090 XI - Xtl) 

4100 XT * X! 2) 

4110 XI2 = XI3) 

4120 Y I = Y ( 1 I 
4130 Yli = Y! 2) 

4140 Y 12 * Y< 3) 

4150 INTRP2 » 0-0 
4160 INTRP4 = 0.0 
4170 INTR2 = 0.0 
4180 INTR4 = 0-0 

4190 CALL PARAI XI, YI, XT« YI 1*XI2« YI 2 »AP1 • AP2 *AP3) 

4200 INTR2 = APl*XT*XT*XT /3.0 + AP2*XT*XT/2. 0 ♦ AP3*XT 
4210 INTR4 = APl*XT**5/5.0 ♦ AP2*XT**4/4. 0 ♦ AP3*XT**3/3.0 
4220 WRITE! 6* 2140) XI II* YI 1) ,STDI 1) 

4230 WRITE! 6* 2210) XI 2) • YI 2 ) . STDI2I • INTR2, INTR4 
4235 IFUD * I POL - 1 
4240 00 4410 I = 3. IFUD 
4250 XP = XI 1-2) 

4260 XT = XII-1) 

4270 XA * XII) 

4280 XAA = XI 1+1) 

4290 YP = YII-2) 

4300 YT = YII-1) 

4310 YA = YII) 

4320 YAA = YII+1) 

4330 CALL PARAI XP, YP.XT, YT ,XA . YA, ABEFl »ABEF2 , ABEF3) 

4340 CALL PARAI XT, YT, XA,YA,XAA •YAA.AAFTl ,AAFT2, AAFT3) 

4350 INTRP2 = 1 1 ABEF1+AAF Tl>/6.0)*l XA**3-XT**3) ♦ 1 1 ABEF2 +AAFT2 )/4.0 1 
1*1 XA*XA— XT*XT )♦ 1 1 ABEF3*AAFT3) /2-0) *1 XA— XT) 

4360 INTRP4 = 1 1 ABEF1+AAFT1 ) / 10. 0) * IXA**5— XT**5) ♦ 1 1 ABEF2+AAFT2 )/0-O)* 
1 1 XA**4— XT**4) ♦ I ( ABEF3+AAF T3) /6« 0) *! XA**3-XT**3 ) 

4370 INTR2 = INTRP2 ♦ INTR2 
4380 INTR4 = INTRP4 ♦ INTR4 
4390 INTRP2 * 0.0 
4400 INTRP4 = 0.0 

4410 WRITEI6.2210) XI I ) , YI I ) ,STD 1 1) • INTR2, INTR4 
4420 XP = XIIPOL-2) 

4430 XT = XtIPOL-1) 

4440 XA = XIIPOL) 

4450 XAA = XI IPOL+1 I 
4460 YP = Y 1 1 POL— 2 ) 

4470 YT = YII POL~ 1 I 
4480 YA = YII POL ) 

4490 YAA * YI IPOL+1) 

4500 CALL PARAIXP.YP, XT, YT,XA,YA, ABEFl, ABEF2.ABEF3) 

4510 CALL PARAI XT *YT »XA»YA,XAA, YAA, AAFT1 ,AAFT2«AAFT3) 
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4520 INTRP2 = ( ( ABEF1 +AAFT l) /6 . 0 ) * ! XA**3-XT**3) + ( { ABEF2 + AAFT2 ) /4.0 ) . 

1*( XA*X A-XT*XT )♦< ( ABEF3+AAF T3) /2.0) * ( XA-XT) 

4530 INTRP4 = !( ABEF1+AAFT l ) / 10. 0) * ( XA**5-XT**5 » ♦ ( ( A8EF2 +AAFT 21/8.0)* 
1! XA**4-XT**4) ♦ { (A8EF3+AAFT3) /6. 0) *( XA**3-XT**3 ) 

4540 INTR2 = INTRP2 + INTR2 
4550 INTR4 = INTRP4 + INTR4 
4560 INTRP2 = 0.0 
4570 INTRP4 * 0.0 

4580 WRITE (6,2250) X( IPQL ) , Y ( I POL ) , STO (I POL) , I NTR2, 1 NTR4 

4590 PARVl = INTR2 

4600 PARV2 = INTR4 

4610 ASSIGN 4630 TO ICON 

4620 GO TO 8000 

4630 INTR2 = INTR2 + INTRP2 

4640 INTR4 = INTR4 + INTRP4 

4650 AGEP = INTR4/(6.0*INTR2) 

4660 WRITE! 6, 2320) 

4670 00 4700 ILLK = IPOL.IDATA 

4680 APR = ALPHA*EXP ( 8ETA*X( ILLK) J 

4690 POEV =( IY( ILLK)-APR)/STD( ILLK) )*100.0 

4700 WRITE (6,2370) X ( ILLK ), Y( ILLK ), STD( ILLK ), APR, POE V 

4710 EXLEN = 1.0/BETA 

4720 WRITE! 6,4730) AGEP, AZERO, EXLEN, INTRP2 ,INTRP4 

4730 FORMAT ( 43X, 36H AGE COMPUTED BY PARABOLIC OPTION IS/ 45X.2PEL5.7 

I, 1 1HS0UARE CMS.///20X, 8H AZERO =, 2PE15.7,20X,20H RELAXATION LENG 
2TH =,2PE15.7//15X,31H INTEGRAL 1 OVER EXPONENTIAL IS,20X, 

331H INTEGRAL 2 OVER EXPONENTIAL I S/22 X ,2PE 15 . 7 ,34X , 2PE15 .7 ) 

4740 IF(0PT7)4750, 4780,4790 
4750 WRITE! 6,4760) 

4760 FORMAT! 18H ERROR CONDITION 6) 

4770 STOP 
4780 GO TO 1012 
4790 0PT7 = 0PT7-1 
4800 GO TO 7000 
5000 FARD = 0.0 
5010 FLUP = 0.0 
5020 FLURG = 0.0 

5030 FARD = X ( IPOL ) /( FLOAT ( I POL- 1 ) ) 

5040 FLUP = FARD*FARD/ 12.0 
5050 CALL POYFZ ( X, Y, IPOL , 2, A ) 

C DIMENSION A! 20) 

5060 FLURG = 2.0 * A(3) * X(IPOL) *FLUP 

5070 ZARD = 0.0 

5080 ZLURG = 0.0 

5090 ZARD * FLUP 

5100 CALL PQYFZ(X*Y2» IPOL, 2, A) 

5110 ZLURG = 2.0* A( 3) * X!IPOL)*ZARD 
5120 ZUCH * FL URG/TRA 1 
5130 ZOOH = ZLURG/TRA2 

5140 ERR0R1 = SORT I ZUCH*ZUCH+ZOOH*ZOOH) 

5150 ERRR = (1.0/6.0)*(1.0/SUM1)*SQRT( ( ZLURG ) **2* (ERRE2)**2*( (6.0*AGET) 
1**2 )*l (FLURG) **2*1 1NTRP2*ERREX )**2) ) /AGE T 
5160 DELAGE = AGET*ERRR*0.5 

5170 WRITE! 6, 5180 ) ZUCH, ZOOH ,ERR0R1 ,DAZ,DBE T , ERREX , ERRR2 , ERRR , DELAGE 
5180 FORMAT! 1H1,54X,9H OPTION 5///l4X,31H FRACTIONAL ERROR IN INTEGRAL 

II, 28X , 3 1H FRACTIONAL ERROR IN INTEGRAL 2/18X.23H OVER TRAPEZOIDAL 


18 


2 RANGE, 37X.23H OVER TRAPEZOIDAL RANGE/23 X.E 13. 7 ,45X, El 3. 7/38X , 

340H FRACTIONAL ERROR IN CURVE I INTEGRATION/ SIX, E13.7///20X, 

417H ERROR IN AZERO = ,E 13.7, 30X, 16H ERROR IN BETA = .E13.7//14X, 

531H FRACTIONAL ERROR IN INTEGRAL 1,28X,31H FRACTIONAL ERROR IN INT 
6EGRAL 2/ 18X. 23H OVER EXPONENTIAL RANGE, 37X, 23H OVER EXPONENTIAL 
7 RANGE/23X, EI3»7,45X,E 13. 7//38X»42H FRACTIONAL ERROR IN THE TOTAL 
8INTEGRATI0N/56X,EI3.7////30X,40H ERROR IN AGE IS EQUAL TO PLUS OR 
9MINUS • 2PE13.5, 12H SQUARE CMS.) 

5190 GO TO 3000 
6000 KORD » ORO + l 

6010 CALL POYFZI X,Y* IPOL* KORD, ARBI) 

6020 MAP * ARB1I KORD+1 I 
6030 LORO = 1 
6040 DO 6050 I = 1,0RD 
6050 LORD = LORD*! KORD- I ) 

6060 FLORD * FLOAT! LORD) 

6070 HEL = X! IPOL) /FLOAT 1 1 POL— I ) 

6080 0RDP2 * ORD+2 
6090 0RDP1 = ORD + 1 
6100 HEL = HEL**0RDP2 
6110 ENY = 1.0 
6120 DO 6130 I = 1,0RDP1 
6130 ENY = ENY*I FLOAT! 1)1 

6140 RERR = HEL*FLORD*X! IPOLl*WAP/( 1 2. 0**ORD) *ENY) 

6150 FSU 1 = RERR/INTR2 

6160 RERR2 * RERR*X! IPOL ) *XI IPOL ) /FLOAT ( « KORD+ll *KORD) 

6170 FSU2 = RERR2/INTR4 

6180 ERFPi = SQRTIFSU1*FSU1 ♦ FSU2*FSU2) 

6190 ERFPR = 1 1.0/6. 0)*! l.O/HSUl )*SQRTURERR2 )**2«- ! ERRE2) **2* I ( 6.0*HAGE 
l)**2f*l ! RERR )**2+! INTRP2*ERREX 1**21 1 /HAGE 
6200 DEHAGE = HAGE*0.5*ERFPR 

6210 WRITE1 6, 62201FSU1 ,FSU2, ERFPI, DAZ.DBET,ERREX,ERRR2. ERFPR, DEHAGE 
6220 FORMAT! 1H1,54X,9H OPTION 6///14X.31H FRACTIONAL ERROR IN INTEGRAL 
11.28X, 31H FRACTIONAL ERROR IN INTEGRAL 2/19X,22H OVER POLYNOMIAL 
2RANGE, 35X.22H OVER POLYNOMIAL RANGE/ 23X,E13.7,45X, 

3E13.7/38X.40H FRACTIONAL ERROR IN CURVE I I NTEGRAT 10N/51X.E13. 7/// 
420X, 17H ERROR IN AZERO =,E 1 3. 7 ,30X , 16H ERROR IN BETA =,E 13.7//14X, 
531H FRACTIONAL ERROR IN INTEGRAL 1 ,28X,31H FRACTIONAL ERROR IN INT 
6EGRAL 2/ 18X,23H OVER EXPONENTIAL RANGE, 37X,23H OVER EXPONENTIAL R 
7ANGE/23X, El 3. 7, 45X, E 13. 7//38X, 42H FRACTIONAL ERROR IN THE TOTAL IN 
8TEGRAT I0N/56X.E 13.7////30X, 40H ERROR IN AGE IS EQUAL TO PLUS OR MI 
9NUS , 2PE13.5, 12H SQUARE CMS.) 

6230 GO TO 4000 
7000 BARD = 0.0 
7010 BLURG * 0.0 

7020 BARD = X ! IPOL )/! FLOAT 1 1 POL- 1) ) 

7030 BLUP = BARD**4/ 180.0 
7040 CALL POYFZ!X,Y, IPOL, 4, A) 

7050 BLURG * BLUP*24.0*A 1 5)*X ( IPOL ) 

7060 CARD = 0.0 

7070 CLURG = 0.0 

7080 CARD = BLUP 

7090 CALL POYFZ t X » Y2» IPOL , 4, A ) 

7100 CLURG = 24.0*A! 5)*CARD*X! IPOL) 

7110 PUCH = BLURG/PARV 1 
7120 POOH = CLURG/PARV2 
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7130 ERR0R3 = SORT l PUCH*PUCH ♦ P00H*P00H1 

7140 ERRRP = ( 1.0/6.0)*! 1 .0/INTR2) *SQRTI ICLURG) **2+ (ERRE2) **2 +( C6.0*AGE 
IP)**2)*1(BLURG)**2*C INTRP2*ERREX) **2) ) /AGEP 
7150 DE AGEP = AGEP*ERRRP*0. 50 

7160 WRITE! 6, 7 1 70 ) PUCH , POOH »ERR0R3 »DAZ ,D BET , ERREX, ERRR2 » ERRRP , DE AGEP 
7170 FORMAT ( 1HL,54X, 9H OPTION 7///14X.31H FRACTIONAL ERROR IN INTEGRAL 
11.28X, 31H FRACTIONAL ERROR IN INTEGRAL 2/15X.30H OVER SMOOTHED PA 
2RAB0LIC RANGE, 31X,30H OVER SMOOTHED PARABOLIC RANGE/2 3X , El 3.7, 45X, 
3E13.7/38X,40H FRACTIONAL ERROR IN CURVE I I NTEGRAT ION/ 51X, E 13. 7/ // 
4 2 OX , 17 H ERROR IN A2ER0 *,E13.7,30X,16H ERROR IN BETA * »E 13 . 7// 14X , 
531H FRACTIONAL ERROR IN INTEGRAL 1,26X,31H FRACTIONAL ERROR IN I NT 
6EGRAL 2/ 18X.23H OVER EXPONENTIAL RANGE ,37X,23H OVER EXPONENTIAL R 
7ANGE/23X, El 3. 7, 45X, E 13. 7//38X, 42H FRACTIONAL ERROR IN THE TOTAL IN 
8TEGRAT I0N/56X,E 13.7////30X, 40H ERROR IN AGE IS EOUAL TO PLUS OR MI 
9NUS , 2PE13. 5, 12H SQUARE CMS.) 

7180 GO TO 1012 

8000 AZERO » 0.0 

8010 BETA * 0.0 

8015 GLQCK = 0.0 

8020 DO 8030 M = IPOL, I0ATA 

8025 WEIGHIM) = 1.0E+04/I STD(M)*STDIM) ) 

8030 GLOCK = GLOCK ♦ WEIGHIM) 

8032 DO 8034 M = IPOL.IDATA 
8034 WEIGHIM) = WE IGHI M J/GLOCK 

8040 CALL WEXPFT(X,Y, WEIGH, IPOL.IDATA, AZERO, BETA, DAZ,D8ETI 
8045 B2 = XI IPOL )*X I IPOL) 

8050 INTRP2=(— AZ ERO/BETA l*EXP( 

1BETA*XI IPOL) ) 

8060 INTRP4=AZER0*I- B2*EXP(BETA* 

1X1 IPOL ) )/BETA + ( 2.0* EXP I BETA* XI I POL) ) / 

2BETA**3)*IBETA*Xl IPOL)- 1.0) ) 

8070 ERREX = SORT! ( DAZ/AZER0)**2+IDBET/8ETA) **2 + I X I IPOL) *DBET ) **2 ) 

8080 DELFN = ( DAZ* INTRP4/AZER0 )**2+ ( XI IPOL ) *1 NTRP4*0BET ) **2+ 1 DBET * AZERO 
1 *EXP( BETA*X ( IPOL) )/I BETA**3)*(2.0*BETA*B2-2.0*Xl IPOL) ) )**2 
8090 ERRE2 = SORT (DELFN) 

8100 ERRR2 = ERRE2/ INTRP4 
8120 GO TO ICOM, (2280, 3640,4630) 

END 
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non* 


IBFTC FTES 


F TEST PROGRAM 

FUNCTION FTESIK,NEST,PE) 

DIMENSION FI 18, 7) 

IFIK-1.EQ.0) GO TO 5 

DATA! IF! I,J), I = 1, 1 8 ) , J=l, 7 ) / 1 . 00, .667 ,.585,.549 ,.528 , .5 15 , .506, .49 

19. . 494.. 490.. 486. .484. .472 ..469. .466.. 463.. 461 . .45 8. . 172 . .133. . 122 

2, .117.. 113,. Ill ..110, .109, .108, 2*. 107,. 106, 2*. 104, 2*. 103, 2*. 102,. 0 

325. . 020. .019.. 018. 7*. 017, 7*. 016,. 0062 ,.0046,. 0044, 2*. 0043, 2*. 0042, 

44*. 0041, 5*.0040, 2*. 0039,. 0015,. 0013,. 0012, 3*. 001 1,9*.0010,3*. 00099 

5. . 00025. .0002. .00019. .00018. 6*. 00017,8*.00016,. 000062, .00005, .0000 

646. . 000044. 2*. 000043, 3*. 000042, .000041, 4*. 000040, 4*. 000039/ 

GO TO 340 

5 IFf NEST. GT. 120) GO TO 360 

IF! PE. GT. 0.995. OR.PE.LT. 0.50) GO TO 380 
10 00 30 I = 1,12 
20 IFINEST.EO.I ) GO TO 160 
30 CONTINUE 

40 IF! NEST .LE. 15) GO TO 100 
50 IFINEST.LE.20) GO TO 110 
60 IF! NEST. LE. 30) GO TO 120 
70 1FINEST .LE.40) GO TO 130 
80 {FINEST. LE.60) GO TO 140 
90 IF! NEST .LE. 120 ) GO TO 150 
100 I = 13 

GO TO 160 
110 I = 14 

GO TO 160 
120 I * 15 

GO TO 160 
130 1-16 

GO TO 160 
140 I = 17 

GO TO 160 
150 I * 18 

160 IF! PE. EO. 0.50) GO TO 250 
170 IFIPE.LE.0.75) GO TO 260 
180 IFIPE.LE.0.90) GO TO 270 
190 IFiPE.LE.0.95) GO TO 280 
200 IF1PE.LE. 0.9751 GO TO 290 
210 1FI PE.LE.0.99) GO TO 300 
220 IF I PE.LE. 0.995 ) GO TO 310 
230 WRITE! 6, 390 ) 

240 STOP 
250 J = 1 

GO TO 320 
260 J * 2 

GO TO 320 
270 J * 3 

GO TO 320 
280 J - 4 

GO TO 320 
290 J * 5 

GO TO 320 
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300 J * 6 

GO TO 320 
310 J = 7 
320 FTES = FU,J) 

330 RETURN 
340 WRITE!6,350) 

GO TO 385 

350 FORMAT! 36H ONLY ONE S1DE0 F TEST NU 1 » 1 ONLY) 

360 WRITE! 6» 370 ) 

GO TO 385 

370 FORMAT! 34H NU 2 IS TOO LARGE NU 2 = 120 MAX.) 

380 WRITE! 6* 390 ) 

385 STOP 

390 FORMAT 1 92H LEVEL OF SIGNIFICANCE OUT OF RANGE 0.50 LESS THAN OR £0 
1UAL TO PE LESS THAN OR EQUAL TO .995) 

ENO 
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. SIBFTC WEXPFT 

C WEIGHTED EXPONENTIAL FIT PROGRAM 

SUBROUTINE WEXPFT ( X, Y, W, M.N, ALPHA. BETA .OAZ.DBET) 

DIMENS IONXl 100 ) » Y I 1 00 ) * W ( 100 ) , H( 1 00 ) 

DOUBLE PRECISION ZUCK .ZACK .ZORP, AKE ,AKK, ARP, FARGO, APPLE 
NAG = 0 
N=N 
M = M 

DO 11 K = M,N 
IF(Y(K))12, 11, 11 

12 YCK) * ABSIY(K)) 

NAG = 1 

11 CONTINUE 

IFINAG) 13*60.13 

13 WRITEI6.6) 

6 FORMAT 1 77HJNEGATIVE Y VALUES HAVE BEEN MADE POSITIVE IN CALCILAIIO 
INS BEYOND THIS LINE. ///I 

60 DO 130 I = M.N 
H(I) = ALOG(YII)) 

130 CONTINUE 

ZUCK * O.ODO 

ZACK = 0.000 

ZORP = O.ODO 

AKE = O.ODO 

AKK = 0.000 

ARP = O.ODO 

DO 10 I - N*N 

ZUCK * ZUCK ♦ Mill 

ZACK a ZACK ♦ X(l)*W(I) 

ZORP * ZORP ♦ H(I)*W(I) 

AKE » AKE ♦ X< I )*W( I I 
AKK « AKK ♦ W(I)*X<I)**2 

10 ARP = ARP ♦ X(I)*HfI)*WU) 

FARGO = ZUCK*AKK - ZACK*AKE 
IF I FARGO ) 30*20*30 

20 WRITE (6*8) 

8 FORMAT ( 30H DIVISION BY ZERO NO SOLUTION.) 

RETURN 

30 APPLE = ( ZORP ♦AKK -ZACK*ARP) /FARGO 
ALPHA * DEXPI APPLE) 

BETA a I ZUCK* ARP -ZORP*AKE) /FARGO 
WRSO a 0.0 

40 DO 50 I * M.N 

50 WRSO » wm*UAPPLE+BETA*X(I) )-H(I ) )**2+ WRSO 

52 WRSO a WRSO /ZUCK 

61 DBET * SORT I ZUCK/FARGQ)*SQRT( WKSQ/FLOAT ( N~M) ) 

70 DAZ a ALPHA* SORT! AKK/FARGO) *SORT ( WRSQ/ FLOAT € N — H ) ) 

RETURN 
END 
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$ I BFTC PARA 

SUBROUTINE PARAiXlt Y1,X2. Y2.X3, Y3» A1,A2 .A3» 

C Y = A 1*X**2 + A2*X *■ A3 

DOUBLE PRECISION XD1,XD2»XD3» YD1 * YD2 # Y03 •DETDl *DETD2 * DETO» OETNl 1 • 
1DETN12* DETN 1* DETN21 » DETN22»D£ TN2» DE TN31 »OETN32 »DET N3 
XD1 = XI 
XD2 = X2 
XD3 = X3 
YD1 = Y 1 
YD2 * Y 2 
YD3 * Y3 

DETDl * XDI*XD1*XD2 + XD1*XD3*XD3 ♦ XD2*XD2*X03 
0ETD2 = XD2*XD3*X03 + XDI*XDl*XD3 ♦ XDi*XD2*XD2 
DETD = DETDl - DETD2 

DETN 1 L = YD I*X02 ♦ YD2*XD3 + YD3*XDl 
DETN12 = XD2* YD3 ♦ XD3*YDI +XD1*YD2 
DETN1 = DETN11 - DETN12 
Al = DETNI/DETO 

DETN2I = XD1*XD1*YD2 ♦ YD1*XD3*XD3 ♦ XD2*X02*YD3 
DETN22 = YD2*XD3*XD3 ♦ YD3*XD1*XDI + YD1*XD2*XD2 
DETN2 * DETN21 - DETN22 
A2 = DETN2/DETD 

DETN31 = XDI*XD1*XD2*YD3 ♦ XD 1+YD2* XD3*XD3 ♦ YD1*X02*XD2*XD3 

DETN32 = YDl*XD2*XD3*XD3 + YD2 *XD 1* X01*XD3 + YD3*XDl*XD2*XD2 

DETN3 = DETN3I - DETN32 

A3 - DETN3/DETD 

RETURN 

END 



• jUBFTC HOMAGE 

C HOW TOE PROGRAM INPUTS AND OUTPUTS. 

SUBROUTINE HOMAGE 

1000 FORMAT ( 1H1, 2X . I 3C 1H* ,3X I , 12HPR0GRAM-ELLI .13 ( 3X ,1H*)///5X , 109HA PRO 
XGRAM FOR THE COMPUTATION OF NEUTRON AGE FROM EXPERIMENTAL ACTIVATI 
XON DATA. SEVEN OPTIONS ARE AVAILABLE.) 

1002 FORMAT I 20X.27H1. TRAPEZOIDAL QUADRATURE ./20X.57H2. POLYNOMIAL Lfc 
X AST— SQUARES FIT AND ANALYTIC QUADRATURE. /20X.51H3. AVERAGED PARAB 
XOLIC FIT AND ANALYTIC QUADRATURE. /2 OX ,97H4. CHOICE OF DIVISION PO 
XINT BETWEEN OPTIONS 1.2. AND 3 QUADRATURES AND WEIGHTED EXPONENT I A 
XL FIT./20X.35H5. ERROR OF OPTION ONE QUADRATURE./ 20 X, 35 H6. ERROR 
X OF OPTION TWO QUADRATURE . / 20X » 37H7. ERROR OF OPTION THREE QUADRA 
XTURE.///I 

1004 FORMAT! 54X» 12HINPUT— FORMAT//5X.61HFIRST CARO BLANK (OTHERWISE THIS 
X CODE DESCRIPTION IS PRINTED) /5X.90HSEC0N0 CARD F0RNATI711) OPTION 
XS 1 THRU 7. 0 IF OPTION SKIPPED, l IF OPTION TO BE EXECUTED.) 

1006 FORMAT (5X.65H0PT ION 1 = TRAPEZOIDAL QUADRATURE AND OPTION 5 ERROR 
XIN OPTION 1 ./5X.76HQPT ION 2 » LEAST-SQUARES ANALYTIC QUADRATURE AN 
XD OPTION 6 ERROR IN OPTION 2./5X.72H0PTI0N 3 * AVERAGED PARABOLIC 
X QUADRATURE AND OPTION 7 ERROR IN OPTION 3./5X.51H0PTI0N 4 = CHOICE 
X OF DIVISION POINT SEE IPOL BELOW.) 

1008 FORMAT! 5X.115HTHIRD CARD F0RMAT(A6.4X,A6,4X,8A6) TYPE OF DATA. UNl 
XTS* ANO LABEL. TYPE OF DATA = AS OR ASR**2 ACCORDING TO WHICH / 
X1X.112H1S READ IN. UNITS * INCH IF R IN INCHES OR SIX BLANK SPAC 
XES IF R IN CMS. LABEL MAY BE ANY ALPHAMERIC SYMBOL. ) 

1010 FORMAT! 5X, l 10HF0URTH CARD FORMAT! I 8.L2X.F 10. 3) NO. OF DATA POINTS 
XANO LEVEL OF SIGNIFICANCE BETWEEN (INCLUSIVE) .50 ANO .995/ 
X5X.92HF IFTH CARD ETC. FORMATOF15.4) R, AS OR ASR»*2, AND SIGMA! AS 
X OR ASR**2) ACCORDING TO CARD 3./ 

X5X. 65HL AST CARD ! IF OPTION 4) FORMAT! 18) INDEX OF DIVISION POINT ! 
X IPOL ) ) 

1012 FORMAT! 5X.85HF0R MORE THAN ONE SET OF DATA REPEAT SEQUENCE STARTIN 
XG WITH CARD TWO OF INPUT-FORMAT.) 

WRITE ( 6* 1000 ) 

WR1TE( 6. 1002 ) 

WRITE! 6. 1004) 

WRITE! 6. 1006) 

WRITE! 6. 1008) 

WRITE! 6. 1010) 

WRITE1 6.1012) 

STOP 

END 


25 


$IBFTC POYFZ 

C POLYNOMIAL LEAST SQUARES FIT 

SUBROUTINE POYFZ (X»Y»N,M,A) 

DIMENSION A( 20) ,Bl 12, 13) .SUMI21), V(12) ,X(100) ,Y( 100) 

COMMON IVDZEO 

IVOZEO = 0 

LS » 2*M ♦ l 

LB = M ♦ 2 

LV * M ♦ l 

DO 5 J = 2,LS 

5 SUMUI = 0.0 
SUM ( 1 ) = N 

DO 6 J = l » L V 

6 VIJ) * 0.0 

00 16 I * UN 
P * 1.0 

vt n = vii) + Yin 
DO 13 J * 2*LV 
P * XI I )*P 
SUM(J) * SUMIJ) ♦ P 
13 VIJ) * VIJ) ♦ YU) * P 
DO 16 J * LBtLS 
P * X I I ) * P 

16 SUMIJ) = SUMIJ) ♦ P 

17 DO 20 I * l.LV 
DO 20 K = 1, LV 
J * K ♦ I 

20 BIK.I) * SUM(J-l) 

DO 22 K = .1, LV 

22 BIK.LB) = VIK) 

23 00 31 L - l.LV 
DIVB * B(L.L) 

IFIDIVB) 25. 42. 25 

25 DO 26 J * L .LB 

26 8IL.J) = B(L.J)/DIVB 
II = L ♦ 1 
IF(Il-LB) 28,33.33 

28 DO 31 I * I l.LV 
FMULTB * BII.L) 

DO 31 J * L , LB 

31 B'( I , J ) * BII.J) - B(L,J)*FMULTB 
33 A I LV ) = BILV.LB) 

1 * LV 

35 SIGMA * 0.0 

00 37 J * I , LV 

37 SIGMA = SIGMA ♦ B( I- 1 » J )*A I J ) 

1 * I- 1 

All) = B I I, LB) - SIGMA 

40 IF I 1-1 ) 41, 41,35 

42 WRITE 16,43) 

224 WRITE (6,225) DIVB,L 

225 FORMAT I E15.7, 16) 

43 FORMAT 1 45H DIVISION BY ZERO. OATA WILL NOT FIT CURVE. 
IVDZEO = 1 

41 RETURN 
END 



SAMPLE PROBLEM AND OUTPUT 


SDATA 

1110111 

AS TRIAL RUN 


30 .05 

.75 


3.81 

21607 . 

1600 . 

5.08 

20174 . 

2600 . 

6.35 

17112 . 

3400 . 

7.62 

13479 . 

4000 . 

8.89 

9970 . 

4000 . 

1 C . 16 

7683 . 

6100 . 

11.43 

549 C • 

5400 . 

12.70 

4003 . 

4850 . 

13.97 

2971 . 

5850 . 

15.24 

215 C . 

7500 . 

16.51 

1525 . 

6250 . 

17.78 

1184 . 

5600 . 

19.05 

887 . 

4800 . 

20.32 

675 . 

8250 . 

21.59 

501 . 

7000 . 

22.86 

375 . 

5200 . 

24.13 

290 . 

2900 . 

25.40 

215 . 

5150 . 

26.67 

178 . 

5650 . 

27.94 

136 . 

6200 . 

29.21 

107 . 

3400 . 

30.48 

90.7 

1850 . 

31.75 

65.3 

3000 . 

33.02 

52.3 

4400 . 

34.29 

34.2 

3300 . 

35.55 

28.7 

10000 . 

36.83 

23.2 

5450 . 

39.37 

13.3 

3600 . 

40.64 

11.2 

5800 . 

43.18 

10.8 

6150 . 


AGE MEASUREMENT BY FOIL ACTIVATION. 

OPTION l LABEL ’TRIAL HUN 


2 


R 

A R 


ST. DEVIATION 

PARTIAL SUM 

PARTIAL SUM 


S 




□F INTEGRAL 

l 

UF INTEGRAL 

2 

0 . 

0 . 


23-2257595E 

03 





38.0999994E-01 

3 1* 3649364E 

04 

23.22 57595E 

03 

39.8334684E 

04 

57. 8226590E 

05 

5.0800000E-00 

52. 06 l 8305E 

04 

67. 096638 7E 

03 

92 .80946456 

04 

17. 204804 2E 

06 

6 3. 499 9 9 95 E-0 l 

68.9998598E 

04 

13. 7096496E 

04 

16.9683635E 

05 

43. 4034758E 

06 

76.1999989E— 01 

78. 2650023E 

04 

23. 2257595E 

04 

26.3196819E 

05 

09.9277468E 

06 

88. 8999996 E- 01 

78. 7950029E 

04 

31.6L28395E 

04 

36 . 29299216 

05 

15. 8328328E 

07 

1.0160000E 01 

79* 3082266E 

04 

62.9676147E 

04 

46. 3325467E 

05 

24. 9857059E 

07 

11. 4299999E 00 

71.7240477E 

04 

7 . 0548244E 

05 

55 .92309 52E 

05 

36. 1 344132E 

07 

12.6999999E 00 

64.5643845E 

04 

78.2256479E 

04 

64.5774097E 

05 

48.697233 7E 

07 

13. 9699999E 00 

57. 9823012 E 

04 

11-4 169 122E 

05 

72.3591242E 

05 

62.495444 8E 

07 

15.2399999E 00 

49.9353828E 

04 

17.4 193 1 96E 

05 

79.21 1 89 69E 

05 

77.04567 53E 

07 

16.5100000E 00 

41, 56 84643E 

04 

17.0362 558E 

05 

85.0223904E 

05 

91 .60534 10E 

07**** 




EXPONENTIAL REGION 

CALCULATIONS 





2 


R 


A R 


ST. DEVIATION 

APPROXIMATION 

PERCENTAGE DEV. 

16.5L0000QE 

00 

4 1. 56 8464 3E 

04 

17.0362558E 

05 

43.43829926 

04 

- 1 .097562 OE— 00 

1 7 . 7800000E 

00 

3 7. 4296 02 IE 

04 

17.7 03190 1 E 

05 

37. 12441066 

04 

17. 239344 8E-02 

19.0500000E 

00 

32. 1894512E 

04 

L7.4L93196E 

05 

31.7282662E 

04 

26.4754746E-02 

2 .03200006 

01 

27.8 7091 16E 

04 

34.0644469E 

05 

27. 1164668E 

04 

22.1475754E-02 

2 1 .58999996 

00 

23.3530173E 

04 

32.6289659E 

05 

23.1 750069E 

04 

54 .55598696— 03 

22 . 8599999E 

00 

19. 5967345E 

04 

27. 1 741 38 5E 

05 

19.Q064494E 

04 

—77.1 744928 E— 03 

24. 1 299999E 

00 

16.8854494E 

04 

16. 8 854494E 

05 

16.92752276 

04 

-24.91674616—03 

25 • 3999999E 

00 

13. 8709396E 

04 

33. 2257390E 

05 

14.46 70558E 

04 

- L 7. 94 1397 7E— 02 

26 . 6699998E 

00 

1 2. 6609422E 

04 

4. 0 187 82 IE 

06 

12 .36422606 

04 

73.8 32344 IE— 03 

27.9399998E 

00 

I.0616753E 

05 

48.3999014E 

05 

1.0567048E 

05 

1.0269524E-02 

29. 20999986 

00 

9L.2949762E 

03 

29.0096185E 

05 

9. 031096 IE 

04 

33.9203219E-03 

3.0479999E 

01 

84.26 3054 BE 

03 

17. 1870618E 

05 

77.1 839895E 

03 

41. 1883 38 8E — 02 

3 1. 75000 COE 

00 

65.8264799E 

03 

3.0241 8 75E 

06 

65 -9650602E 

03 

-45.8239036E-04 

33 . 0200000E 

00 

5 7.0237560E 

03 

47.9740963E 

05 

56 • 3768435E 

03 

13.48462256-03 

34.29000006 

00 

4.0212499E 

04 

38.8015342E 

05 

48 • 1822968E 

03 

-2.0539902E-01 

35 . 54999 9 7E 

00 

36. 271 l 301 E 

03 

12.6 380247E 

06 

41.2298169E 

03 

— 39.2362456E— 03 

36.8299999E 

00 

31.4696138E 

03 

73.926463 IE 

05 

35. 193385 IE 

03 

-5.0371291E— 02 

39.3699999E 

00 

2.06 1 4958E 

04 

55 . 7998867E 

05 

25 . 7060034E 

03 

-91.23754696-03 

4 . 0640000E 

01 

1S.4980271E 

03 

95.7933540E 

05 

21.9695585E 

03 

-36.23979386-03 

43 • 1 799998E 

00 

2.0136733E 

04 

1 1.4667509E 

06 

16 .0 4703 78E 

03 

35. 6656866E— 03 


AGE COMPUTED BY TRAPEZOIDAL OPTION IS 
45.3679B05E OOSQUARE CHS. 


AZERO = 33.4691639E 05 


RELAXATION LENGTH -B.O857326E-00 


INTEGRAL 1 OVER EXPONENTIAL IS 
35. 12304746 05 


INTEGRAL 2 OVER EXPONENTIAL IS 
23*54400 L6E 08 


OPTION 5 


FRACTIONAL ERROR IN INTEGRAL 1 
OVER TRAPEZOIDAL RANGE 
0 • 60 19756E— C2 

FRACTIONAL ERROR IN CURVE I 
0.6028 l 50 E- 02 


FRACTIONAL ERROR IN INTEGRAL 2 
OVER f RAPEZO I DAL RANGE 
0. 31800926-03 

INTEGRATION 


ERRDR IN AZERO =0. 17811168 06 


ERROR IN BETA *0. 21 67 722E-Q2 


FRACTIONAL ERROR IN INTEGRAL 1 
OVER EXPONENTIAL RANGE 
0.664 8373E— 0 1 


FRACTIONAL ERROR IN INTEGRAL 2 
OVER EXPONENTIAL RANGE 
0 • 67556 16E-0 i 


FRACTIONAL ERROR IN THE TOTAL INTEGRATION 
0 . 52546 53E-01 


ERROR IN AGE IS EQUAL TO PLUS OR MINUS 11.9I965E-01 SQUARE CMS. 





AGE 

MEASUREMENT 8Y FOIL 

ACT XV AT ION. 





OPTION 2 


LABEL -TRIAL RUN 




R 


2 

A R 

s 


STANDARO DEVIATION 

APPROXIMATION 

0 . 


0 . 


0.2322576E 

05 

-0.6331250E 

02 

0. 38100 OOE 

01 

0. 3136494E 

06 

0.2322576E 

05 

0.3139U8E 

06 

0.5080000E 

01 

0.520 61 83E 

06 

0.&709664E 

05 

0. 5234440E 

06 

0. 6350000E 

01 

0.68999868 

06 

0.1370965E 

06 

C-6943553E 

06 

0.7620000E 

01 

0. 7 8265006 

06 

0.2322576E 

06 

0.77B1188E 

06 

0.889OOCOE 

01 

0.78795 OOE 

06 

0.3161284E 

06 

0. 8053607E 

06 

0.1 Cl 60008 

02 

0. 79 30 82 3E 

06 

0.6296761E 

06 

0. 7 T94863E 

06 

O.U430COE 

02 

0.7172405E 

06 

0. 7054824E 

06 

Q.72Q3U7E 

06 

0.12 700 OOE 

02 

0. 64 5643 BE 

06 

0.7B22565E 

06 

0- 64769096 

06 

0. 1397000E 

02 

0. 579B230E 

06 

0* U41691E 

07 

0.575L425E 

06 

0.15240 OOE 

02 

0.4993538E 

06 

0.17419328 

07 

0.5034831E 

06 

0. 16510006 

02 

0. 4156B46E 

06 

0. i 703626E 

07 

0.414452 7E 

06 


SIGMA 1 SQUARED 


0. 10Q4194E 09 


LEVEL OF SIGNIFICANCE -0.750QE 00 
EXPONENTIAL REGION calculations 


R 


A R 

5 


ST. DEVIATION 

APPROXIMATION 

PERCENTAGE DEV. 

I6.51000COE 

00 

41.56 8 4 643E 

04 

17.0362558E 

05 

43.4382992E 

04 

- 1.09 756 20E-00 

17.7600000E 

00 

37.4296C21E 

04 

17.70319016 

05 

37. 12*41068 

04 

17.23 9344 BE-02 

19.0500000E 

00 

32.1894512E 

04 

17.4193196E 

05 

3 1 • 72B2662E 

04 

26.47547466-02 

2.0320000E 

01 

27.870911 6E 

04 

34.0644469E 

05 

27.11646686 

04 

22.14 75 754E— 02 

21.5899999E 

00 

23.35301 73E 

04 

32.6289659E 

05 

23.1750069E 

04 

54.5559869E-03 

22 . 8599999E 

00 

19. 5967345E 

04 

27. 1 741385E 

05 

19. 80644 94E 

04 

— 77. 1 744928E— 03 

24.1299999E 

00 

16. 88544 94E 

04 

16.8854494E 

05 

16.9275227E 

04 

-24. 9 16746 IE —03 

25. 3999999E 

00 

13-8709396E 

04 

33.22573906 

05 

14.46 70558E 

04 

-17.94139776-02 

26.6699998E 

00 

12.6609422E 

04 

4.01878216 

06 

12.3642260E 

04 

73.83234416-03 

27.93999988 

00 

1.0616753E 

05 

48.3999014E 

05 

1. 056704 BE 

05 

1 • 0269524E— 02 

29.2099998E 

00 

91.2949762E 

03 

29.00961858 

05 

9.0 31096 IE 

04 

33. 92032 19E-03 

3.0479999E 

01 

84.26305488 

03 

17.18 7061 8£ 

05 

77. 1839895E 

03 

41.18833886-02 

31.75000CCE 

00 

65.8264799E 

03 

3.02418756 

06 

6S.9650602E 

03 

—45. #2390368—04 

33-0200000E 

00 

57.0237560E 

03 

47.9740963E 

05 

56.3768435E 

03 

13.4 846 22 5£— 03 

34.2900000E 

00 

4.0212499E 

04 

38.80153426 

05 

48. 1822968E 

03 

-2.05399O2E-0I 

35.5499997E 

00 

36.271 1301E 

03 

12.6380247E 

06 

41*22981696 

03 

— 39.2362456C— 03 

36.8299999E 

00 

31.46961 3 BE 

03 

73.92646318 

05 

35. 1933851E 

03 

—5. 037 129 IE— 02 

39.3699999E 

oo 

2.061495BE 

04 

55.79988678 

05 

25. 70600 34E 

03 

—91 .237546 9€— 03 

A. 06 400 OOE 

01 

18.49802718 

03 

95.79335406 

05 

21.96955858 

03 

-36.23979388-03 

43 . 1 799998E 

00 

2.0136733E 

04 

1 1 .46675098 

06 

16.0470378E 

03 

3 5 .66568668—03 


45.8696127E OOSQUARE CMS. 


A2ER0 - 33.46916396 05 


RELAXATION LENGTH - -6. 0857326E- 00 


INTEGRAL I OVER EXPONENTIAL IS 
35. 12304748 05 


INTEGRAL 2 OVER EXPONENTIAL IS 
23.54400166 08 
DEGREE OF EQUATION * 5 


OPTION 6 


FRACTIONAL ERROR IN INTEGRAL l FRACTIONAL ERROR IN INTEGRAL 

OVER POLYNOMIAL RANGE OVER POLYNOMIAL RANGE 

0. 36 12392E— 06 0. 21432686-07 

FRACTIONAL ERROR IN CURVE I INTEGRATION 
0.361B745E-O6 


2 


ERROR IN A2ER0 -0.178U16E 06 


ERROR IN BETA -0.2I67722E-02 


FRACTIONAL ERROR IN INTEGRAL I 
OVER EXPONENTIAL RANGE 
0.6648373E— 01 


FRACTIONAL ERROR IN INTEGRAL 2 
UVEA EXPONENTIAL RANGE 
0. 67556168-01 


FRACTIONAL ERROR IN THE TOTAL INTEGRATION 
0.5246*186-01 


ERROR IN AGE IS EQUAL TO PLUS OR MINUS 12.03256E-01 SQUARE CMS. 


PERCENTAGE 
DEVIATION 
2 7. 2 5 960 3 5E— 02 
11 -29B4059E— 01 
42.11 352 16E-Q1 
41,16323958-01 
19.50925148-01 
55. 07453 5 8£- 01 
21. 59206 3 9E— 01 
43.5342598E-02 
26. 16 818696-02 
4.0996335E-01 
23* 70527 03E— 02 
72.31 01673E-03 


29 


age measurement by foil activation. 

OPTION 3 LABEL = TRl AL RUN 


2 

R A R 

S 

0 . 0 . 

3 8. 0 99 99 <34 E— 01 3U3649364E 04 

5.0800000E-00 52.0618305E 04 

63.4999995E-01 68.9998598E 04 

76. 1999989E-Q 1 78.2650023E 04 

88.8999996 E— 0 1 78.7950029E 04 

1 .0 160000E 01 79.3082266E 04 

11.4299999E 00 71.7240477E 04 

12.6999999E 00 64.5643845E 04 

13.9699999E 00 57.9823012E 04 

15.2399999E 00 49.9353828E 04 

**** 16.5100000E 00 4 1 . 56 84643E 04 


ST. OE V I AT I ON PARTIAL SUM 

OF INTEGRAL 1 

23. 22 57595E 03 

23. 2257595E 03 45.1170650E 04 

67 . 0966 387E 03 98.0209846E 04 

13.7096496E 04 17.550008IE 05 

23.2257595E 04 26.9881523E 05 

3 l .6 128395E 04 37.0077753E 05 

62.9676147E 04 47.0902677E 05 

7 . 0548244E 05 56.7214184E 05 

78. 2256479E 04 65.3704309E 05 

11. 4 16912 2E 05 73.1568403E 05 

17.4193196E 05 8.0019057E 06 

17.0362558E 05 85.8088722E 05 


PARTIAL SUM 
OF INTEGRAL 2 

36.9945 I21E 05 
14.47067 1 5E 06 
4.01 403 2 4E 07 
86.489601 IE 06 
15. 49 1 4529E 07 
24 . 6536996E 07 
35. 8679900E 07 
48.4462638E 07 
62.2 790976E 07 
76.8942471E 07 
9 1. 45757 01E 07**** 


EXPONENTIAL REGION CALCULATIONS 


2 


R 


A R 

c 


ST. DEVIATION 

APPROXIMATION 

PERCENTAGE DEV. 

16. 5100000E 

00 

o 

41. 5684643E 

04 

17.0362558E 

05 

43.43S2992E 

04 

-1. 09756206-00 

17.78000006 

00 

37.429602 IE 

04 

17.7 031901 E 

05 

37. 1244106E 

04 

17. 239344 8E-02 

19. 0500000E 

00 

32- 18945 12E 

04 

17.4193196E 

05 

31.728 2662E 

04 

26. 4754 746 E— 02 

2.0320000E 

01 

27.8709L16E 

04 

34. 0644469E 

05 

27. 1164668E 

04 

22.14757546-02 

21.5899999E 

00 

23. 35301 73E 

04 

32.6209659E 

05 

23. 1750069E 

04 

54.55598696-03 

22.8599999E 

00 

19. 5967345E 

04 

27.17413856 

05 

19. 8064494E 

04 

— 77.1744928E— 03 

24 . 1299999E 

00 

16. 8854494E 

04 

16. 8854494E 

05 

16.9275227E 

04 

—24.916746 IE— 03 

25.3999999E 

00 

1 3. 8709396E 

04 

33.2257390E 

05 

14.46 705586 

04 

— 17.9413977E-02 

26 . 6699998E 

00 

1 2 . 66 09422E 

04 

4. 0187 82 IE 

06 

12.3642260E 

04 

73.8 32344 IE— 03 

27.9399998E 

00 

1.06167536 

05 

48.3999014E 

05 

l . 056704 8E 

05 

1.0269524E-02 

29.2099998E 

00 

91-2949762E 

03 

29.0096185E 

05 

9.0 3 1096 LE 

04 

33.92032 19E— 03 

3 . 0479999E 

01 

84.263054 8E 

03 

17.18 70618E 

05 

77 . 1039895E 

03 

41.18 83388E— 02 

31. 7500000E 

00 

65. 8264799E 

03 

3.0241875E 

06 

65.96506026 

03 

— 45 .8239036E— 04 

33.0200000E 

00 

57.0237560E 

03 

47.9740963E 

05 

56. 3768435E 

03 

l 3. 4846 225 E— 03 

34.2900000E 

00 

4.02 12499E 

04 

38. 80 15 34 2E 

05 

48. 1822968E 

03 

- 2.05 3 9902 E— 01 

35.5499997E 

00 

36. 271 130 IE 

03 

12.63002476 

06 

41 .2298 16 9E 

03 

— 39.2 3 62 456 E— 03 

36. 8299999 E 

00 

31.46961386 

03 

73. 926463 IE 

05 

35. 1933851E 

03 

- 5. 037 129 IE— 02 

39. 3699999E 

00 

2.06 14958E 

04 

55.7998867E 

05 

25* 7060034E 

03 

-91.2375469E-03 

4 . 0640000E 

01 

1 8. 4980271E 

03 

95.7933540E 

05 

21.9695585E 

03 

— 36.2 39793 8E — 03 

43.1799998E 

00 

2. 01 36 733E 

04 

11.4667509E 

06 

16.0470378E 

03 

35.6656866E— 03 


AGE COMPUTED BY PARABOLIC OPTION IS 
45. 052 563 7E OOSQUARE CMS. 


AZERO = 33. 46916396 05 


RELAXATION LENGTH = -8. 08 57326E- 00 


INTEGRAL 1 OVER EXPONENTIAL IS 
35. 123Q474E 05 


INTEGRAL 2 OVER EXPONENTIAL IS 
23.5440016E 08 


FRACTIONAL ERROR IN INTEGRAL 1 
OVER SMOOTHED PARABOLIC RANGE 
0*16070896-03 

FRACTIONAL ERROR IN CURVE 
0. 1610369E— 03 


FRACTIONAL ERROR IN INTEGRAL 2 
OVER SMOOTHED PARABOLIC RANGE 
0. 1027 274E- 04 

INTEGRATION 


ERROR IN AZERO *0.17811166 06 


ERROR IN BETA *0-21677228-02 


FRACTIONAL ERROR IN INTEGRAL 1 
OVER EXPONENTIAL RANGE 
0. 6648373E— 01 


FRACTIONAL ERROR IN INTEGRAL 2 
OVER EXPONENTIAL RANGE 
0.6755616E-01 


FRACTIONAL ERROR IN THE TOTAL INTEGRATION 
0.52347276-01 


ERROR IN AGE IS EQUAL TO PLUS OR MINUS 11.79189E-01 SQUARE CMS. 
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